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Abstract. We reanalyse optical gravitational lens sur- 
veys from the literature in order to determine relative 
probabilities in the Ao-rio plane, using a softened singular 
isothermal sphere lens model. In addition, we examine a 
portion of the Ao-rio plane which includes all viable cos- 
mological models; this is vital for comparison with other 
cosmological tests. The results are, within the errors, con- 
sistent with those of more specialised analyses, such as 
those concerning upper limits on Aq in a flat universe. 
We note that gravitational lensing statistics can provide 
a quite robust lower limit on the cosmological constant 
as well, which could prove important in confirming cur- 
rent claims of a positive cosmological constant. At 95% 
confidence, our lower and upper limits on Aq — f^Oj using 
lens statistics information alone, are respectively —3.17 
and 0.3. For a flat universe, these correspond to lower and 
upper limits on Aq of respectively —1.09 and 0.65. 
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1. Introduction 

The use of gravitational lensing statistics as a cosmological 
tool was first considered in detail by Turner et al. (1984); 
the influence of the cosmological constant was investigated 
thoroughly by Fukugita et al. (1992), building on the work 
of Turner (1990) and Fukugita et al. (1990). More recently, 
Kochanek (1996, hereafter K96, and references therein) 
and Falco et al. (1998) have laid the groundwork for using 
gravitational lensing statistics for the detailed analysis of 
extragalactic surveys. However, these analyses either have 
concentrated on a small subset of the possible cosmologi- 
cal models as described by the density parameter Qq and 
the cosmological constant Aq, have used a simpler (sin- 
gular) lens model or both. This analysis is the first time 
Ao and JIq have been used as independent parameters in 
conjunction with a non-singular lens model in an analysis 



of this type, complementing similar analyses with other 
emphases. (See Cheng & Krauss (1999) for a discussion of 
the importance of including a core radius.) Also, we in- 
clude enough of the Ao-rio plane to avoid neglecting any 
possibly viable models; this also makes the comparison 
with a variety of other cosmological tests easier. This is 
especially important in light of the fact that many analy- 
ses (e.g. Perlmutter et al. 1998; Riess et al. 1998; Schmidt 
et al. 1998; Carlberg et al. 1998a; Lineweaver 1998; Guerra 
et al. 1998; Daly et al. 1998) are now suggesting that our 
universe may contain a significant cosmological constant 
and be non-fiat. 

The plan of this paper is as follows. Sect. 2 reviews the 
groundwork and serves to define our notation. In Sect. 3 
we specify the observational data and selection functions 
we use and formulate prior information about the parame- 
ters Ao and flo- Sect. 4 describes the parametric submod- 
els we use and the numerical computations we perform. 
In Sect. 5 we discuss our results and compare them with 
others. Sect. 6 presents our summary and conclusions. 

2. Probability of multiply imaged sources 

In this section we briefiy review the statistical concepts 
introduced in K96; this also serves to define our notation. 
Note that with regard to cosmogical notation we follow 
that of Kayser et al. (1997), repeating here only 2 equa- 
tions needed for discussion in this paper: the comoving 
spherical volume element at redshift z reads 



where 



dz 



Q{z) = Qail + zf - {Qo + Ao - 1)(1 + zf + Aq. 



(1) 



(2) 
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Following the K96 approach, we assume that the light 
defiection properties of the gravitational lenses can be 
modelled with a particular type of circularly symmetric 
lens models with a monotonically declining radial mass 
profile. Such lens models generally create three images and 
have two critical radii on which the magnification diverges 
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(e.g. Schneider ct al. 1992). It is possible to estimate the 
probability p{m, Zg) of the event 

A source at redshift Zs is triply imaged. The total 
apparent magnitude of the three images is m. The 
image configuration meets the selection criteria S 
and, particularly, shows the properties C. 

If the outer and inner critical angular radii of the lens po- 
tential are respectively r-|_ and r_ , the image magnification 
at radial angular position r is /i(r), the total magnifica- 
tion of the three images of a source at angular position y is 
M(y), the hmctions S{y) and C{y) are 0|1 valued selection 
functions, the comoving density of lenses of luminosity L 
is dr?,/dL and the number-magnitude counts of sources are 
dVV/dm, then 

r_ 

X B{m,z„y)Siy)C{y)drdLdz, (3) 



where 



dN 



B{m,z,y) = — {m-h2.51og[M(y)],^;} 
dm 



dN . 



•(4) 



The critical radii, the image magnifications and the source 

position are functions of the lens model, the luminosity of 
the lens galaxy and the redshifts of the source and the lens 
galaxy. If the underbraced functions are dropped, Eq. (3) 
yields the optical depth - the fraction of the sky included 
within the caustics of all lenses between us and the sources 
at redshift ^s- The inclusion of these functions accounts 
for magnification bias, survey selection effects (including 
what is defined as a lensing event) and allows the observed 
image separation to be taken into account. 

Equation (3) parametrically depends on Xq and fio 
through Eq. (1) and through the angular size distances,^ 
which are needed for calculating observable quantities 
from the lens model (these also depend on the source 
and lens redshifts). Equation (3) additionally depends on 
parametric submodels required to model the lens popula- 
tion and the number-magnitude counts of sources. Since 
throughout this paper we are principally interested in Aq 
and Oo, hereafter we refer to the submodel parameters 
as nuisance parameters (although technically they are on 
the same footing with Aq and flo, there are not of as 
much interest here and thus a nuisance). In principle, one 
could also incorporate other observables into the paramet- 
ric model; the reasons for not doing so are practical. 



^ In general, the angular size distances depend not only on 
Ao and Qq but on the degree of homogeneity in the universe as 
well (see, e.g., Kayser et al. 1997). However, in contrast to some 
other cosmological tests, this effect is relatively unimportant 
for the type of analysis performed here (see, e.g., Fukugita 
et al. 1992). 



Assuming the siirvey selection function S is known, we 
can numerically compute Eq. (3) and reasonably estimate 
the probability 1 —p{mi,Zi) that the quasar i is singly 
imaged or the probability p(mi, Zi,9i) that the quasar i is 
multiply imaged and its images (within some tolerance) 
are separated by 6i. If the survey data D contains M 
singly and N multiply imaged quasars, we can estimate 
the probability of the event 

In a model universe fixed by the cosmological pa- 
rameters Xq, flo and the nuisance parameters ^, a 
multiply imaged quasar survey collects the observa- 
tional data D. 

by applying the parametric model (or likelihood function) 

M 

ln[p(£)|Oo,Ao,^)] = -'^p{mi,Zi) 



N 



+ ^^^\p{mj,Zj,9j)], 



(5) 



where the logarithm ln[l — p{mi,Zi)] was expanded to 
first order. We can combine surveys of different objects by 
adding the logarithms of the likelihood functions for the 
individual surveys, and can combine surveys containing 
the same objects by applying their joint selection func- 
tion. 

In Bayesian theory the model parameters Aq, Oq, ^ 
are regarded as random quantities with known joint prior 
probability density function p(Ao, Oq, ^). Applying Bayes's 
theorem, the appropriate posterior probability distribu- 
tion given the observational data D is 



p{Xo, no^\D) = p{D\Xo, flo, (3 p(Ao, Oo, ^), 



(6) 



where the operation '(g)' denotes multiplication followed 
by normalisation. Marginalising the nuisance parameters 



p{Xo,flo\D) = J p{Xo,flo,^\D)d^. 



(7) 



yields the (marginal) posterior probability density func- 
tion for the parameters Aq and Oq. In the limit where 
all nuisance parameters take a precise value, ^ = ^o, the 
joint prior probabihty density function p(Ao,Oo,^) fac- 
torises into p(Ao,Oo) and a delta distribution — ^o), 
and Eq. (7) simplifies to 



p{Xo, Oo|r>) = p{D\Xq, Oo, ^o) O p(Ao, Oo). 



(8) 



On the basis of Eq. (7) or Eq. (8), we can calculate confi- 
dence regions for two parameters or perform further mar- 

ginalisations and calculate mean values, standard devia- 
tions and marginal confidence intervals for one parameter. 
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Table 1. Observational data of multiply imaged quasars con- 
tained in the sample. The magnitudes are V magnitudes unless 
otherwise specified. The image separations are taken from Ko- 

chanok ot al. (1997) 



Identifier 


m [mag] 


Zb 


en 


Q 0142-100 


17.0 


2.72 


2.2 


Q 1009-0252 


18.1 B 


2.74 


1.5 


Q 1115+080 


16.2 


1.72 


2.2 


Q 1208+1011 


17.9 


3.80 


0.48 


Q 1413+117 


17.0 


2.55 


1.2 



3. Observational data and prior information 

We use the observational data of the optical multiply im- 
aged quasar surveys by Crampton et al. (1992), Jaunsen 
et al. (1995), Kochanek et al. (1995), Yee et al. (1993) 
and the observational data of the HST Snapshot Survey 
compiled by Maoz et al. (1993), including Q 0142-100, 
Q 1115+080 and Q 1413+117. If applicable, we replace 
the apparent quasar V magnitude catalog data found in 
Crampton ct al. (1992), Jaunsen et al. (1995) and Ycc 
et al. (1993) with more current data from Veron-Cetty 
& Veron (1996). We estimate the Kochanek et al. (1995) 
apparent quasar V magnitude data by adding the survey 
average V-R and V-I colours to the observational R and 
I magnitude data. Following K96, we only include quasars 
with redshift Zs > 1 . In all, our sample contains 807 singly 
and 5 multiply imaged quasars. The observational data of 
the multiply imaged quasars are summarised in Table 1. 
Our complete input data can be obtained from 

http : //multivac . jb . man .ac.uk: 8000/ ceres 
/data_f rom_papers/lower_liinit . html 

This follows K96 for purposes of comparison. Since much 

larger surveys (i.e. CLASS) will be considered in a future 
paper, there is little point in increasing the number of 
lenses for its own sake. Since radio observations are con- 
sidered in more detail in a companion paper (Helbig et al. 
1999), we restrict ourselves to optical surveys in this pa- 
per. We use the Crampton et al. (1992), HST Snapshot 
Survey and Yee et al. (1993) survey selection functions 
proposed in Kochanek (1993), the Jaunsen et al. (1995) 
survey selection function at I'.'O seeing and the prelimi- 
nary Kochanek ct al. (1995) survey selection function. 

Before considering prior information in more detail, 
one must first decide which region of the Aq-Oo plane is 
to be investigated. Clearly, this region should be defined 
by either exact constraints or conservative estimates, as 
opposed to current 'best fit' values (and their errors), in 
order to avoid excluding any possibly viable cosmological 
models. Also, it is desirable for the region to be on the 
large side, so that in addition the sensitivity of the test 
(i.e. what regions of the Ao-$lo plane can be ruled out at 
a high confidence level) can be investigated. 



3 

3.1. The range ofQo and Aq 

The mass clustered with galaxies on smaller scales, Slo.gai, 
is 0.1 within a factor of two (e.g. Peebles 1993). This lower 
limit is small compared to our full flo range so we do not 
assume any prior lower limit on Qq except, of course, 

fio > 0. (9) 

Especially for comparison with other work it is important 
to note that, within the framework of c;osmological mod- 
els based on general relativity with which we (and almost 
everyone else at present) are working, > is a require- 
ment. Results reported which include fio < within the 
errors, or even as a best-fit value, do not indicate 'im- 
plausible results' but merely improper statistics. Often, 
confidence contours are assumed to be ellipses and these 
are extended, if applicable, to < 0. (Of course, it is 
possible that Oq = is within the errors or even the best 
fit value for a certain set of results.) 

An extremely conservative upper limit comes from dy- 
namical tests on larger (though still cosmologically small) 
scales; when this work was started, we assumed an (again, 
extremely conservative) upper limit ^Iq < 2 (Czoske 1995). 
Since then, these methods have started to indicate smaller 
values of ilg, (e.g. da Costa et al. 1998) more in line with 
both a long tradition of low Qq values (e.g. Gott et al. 
1974; Coles & ElUs 1994, 1997) (albeit with somewhat 
larger errors) as well as new determinations (often with 
quite small errors), examples of which are mentioned in 
Sect. 3.2. 

We have assumed no prior upper or lower limits on Aq 
per se. This has two reasons: 

— 'Direct' measurements of Xq (as opposed to measure- 
ments of a combination of parameters involving Aq) 
are virtually nonexistent. 

— We obtain a small enough range in Aq from the values 
obtained from joint constraints on the range of and 
Ao. 

Historically, positive Aq values have been considered 
more than negative ones, probably because positive val- 
ues can have a wide range of relatively easily observable 
effects, while negative ones are more difficult to measure. 
Many cosmological tests have a degeneracy such that Aq 
and fio are correlated, so that increasing Aq can be com- 
pensated for in some sense by increasing flo as well. Thus, 
effects of negative values of Aq for a given value of flo are 
hard to differentiate from the effects of larger values of CIq 
for larger (less negative) values of Aq or even Aq = 0. 

Here, we consider negative values of Aq as well. There 
is no a priori reason why they cannot exist. If one be- 
lieves that the 'source' of Aq are zero-point fluctuations of 
a quantum vacuum, this would lend support to the idea 
that Aq > 0. However, it is not clear that this mMst be the 
only source of Aq, and indeed it has been argued that, if 
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this source of Xq exists, there must be an additional con- 
tribution with a negative value (e.g. Martel et al. 1998, 
though the assumption that this is possible is so obvious 
to the authors it is barely stated!). 

In spatially closed (A; = +1) models, the antipode is re- 
quired to be at 2; > 4.5, the redshift of the most redshifted 
multiply imaged object currently known (Gott ct al. 1989; 
Park & Gott 1997).^ The light grey shaded area in the 
panel in the middle of the left column of Fig. 1 marks 
the right side of the region thus enclosed. This gives us a 
slightly Qo-dependent upper limit on Aq which is slightly 
stronger than that obtained by merely excluding models 
with no big bang. (This can be done because these models 
have a maximum redshift which is less than the redshift of 
high-redshift objects, the only exception being some cos- 
mological models which have fig < 0.05, the robust lower 
limit discussed above (e.g. Feige 1992).) 

The age of the universe in units of the Hubble time. 



To 



IS 



00 

/ 



dz 



{l + z),/Q{z)' 



(10) 



where Q{z) is given by Eq. (2) and thus depends on 
and Aq. (There are world models in which the maximum 
redshift is not infinite but these are all models without 
a big bang and are excluded by the constraint from the 
antipodal redshift or the lower limit on Qq as discussed 
above and are thus not relevant for this work.) Clearly, in 
any physically realistic world model, tqHq^ exceeds the 
age of the oldest galactic globular clusters: 



T0> IscHq. 



(11) 



Following Carroll et al. (1992), we take a robust lower 
limit on Aq from conservative lower limits on the Hubble 
constant and age of the universe. This gives a lower limit 
on Aq from the value at fio = 0; at larger values of Q.q 
the constraint on Aq is not as strict — by assuming the 
lower limit of Aq = — 5 independent of i7o we are being 
conservative. We choose Aq > —5 instead of Aq > —7 as in 
Carroll et al. (1992) since no published current constraints 
examine this region in detail. (Were this the case, then 
including this area would be helpful if only to aid a direct 
comparison.) This value corresponds roughly to the one- 
sided 99% confidence level in the top row of Fig. 1 (see 
Sect. 3.2), which is also a reason not to extend the area 
to more negative Aq values. 



^ Recently, a lensed object of even larger redshift hcis been 
detected z = 4.92 (Pranx et al. 1998). However, at our reso- 
lution this would make only a negligible difference to the results 
so we have not updated the calculations to reflect this. 



3.2. Prior probability for Aq and flo 

We have assumed no prior knowledge of Aq per se, apart 
from the upper and lower limits discussed above. This has 
three reasons: 

— 'Direct' measurements of Aq (as opposed to measure- 
ments of a combination of parameters involving Aq) 
are virtually nonexistent. 

— Based on general knowledge from the literature and 
our own low-resolution calculations, we expect lens 
statistics itself to constrain Aq quite well. 

— Although recent measurements are encouraging (see 
Sect. 5), the value of Aq is observationally not as well 
established as that of fio- 

Regarding tgc and Hq as independent random quanti- 
ties with known prior probability density functions p{tgc) 
and p{Ho), the probability that Eq. (11) is satisfied is 



TO /Ho 



P{tq > tgcHo 



I AH.) I 



p(tgc)digc dHo. 



(12) 



A cosmological world model is compatible with the ab- 
solute age of the oldest galactic globular clusters as long 
as the above expression does not vanish. Reasonably, we 
assume a prior probability density function that is pro- 
portional to this expression 



TO /Ho 



Pi(Ao,Oo) = l«) j p{Ho) j 



p{tgc)dtgc dHo. 



(13) 



The best estimate of the absolute age of the oldest galactic 
globular clusters currently is tgc = 11.5 ± 1.3 Gyr (Cha- 
boyer et al. 1998). We choose to formulate this prior infor- 
mation in the form of a lognormal distribution that meets 

these statistics 

p(tgc) = L(V|ll-5 Gyr, 1.3 Gyr). (14) 

Similarly, we roughly estimate Hq = 65± 10 kms~^ Mpc~^ 
and choose to formulate this prior information in form of 
a normal distribution 

p{Ho) = iV(Fo|65kms-^Mpc-\lOkms-^Mpc-^), (15) 

where the notation for L and N is such that the two ar- 
guments correspond to the mean and standard deviation. 

This estimate is compatible with 'small' values of the 
Hubble constant, which is conservative in the sense that 
it restricts our region of the Ao-Oo plane less than would 
'large' values. By the same token we neglect any time be- 
tween the big bang and the formation of the oldest globu- 
lar clusters. Inserting Eq. (14) and Eq. (15) in Eq. (13) one 
obtains a well-founded a priori probability distribution for 
the parameters fio und Aq. 

Although observational evidence has always indicated 
a low value of VLq (e.g. Gott et al. 1974; Coles & Ellis 1994, 
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1997), the inflationary paradigm (e.g. Guth 1981), coupled 
with a prejudice against a non- negligible value of Aq, has 
created a prejudice in favour of Oq = 1)^ unfortunately 
too often to the extent where this prior belief has been 
elevated to the status of dogma (see, e.g., Matravers et al. 
1995, for an illuminating account) even though there are 
serious fundamental problems with the inflationary idea 
(e.g. Penrose 1989) and even though there might be other 
solutions to the problems it claims to solve (e.g. Barrow 
1995; Collins 1997). What is more, some current infla- 
tionary thinking (e.g. Turok & Hawking 1998) seems able 
to predict values for Aq and CIq similar to current obser- 
vationally determined values, though it would have been 
more interesting had this prediction been made before the 
recent improvements in the observational situation. (To 
be fair, many leading practitioners of inflation consider 
a flat universe to be a robust prediction and its obser- 
vational falsification essentially a falsification of the entire 
paradigm.) Recently, in the light of overwhelming observa- 
tional evidence in favour of a low value of flo (e.g. Carlberg 
et al. 1998b; Carlberg 1998; Carlberg et al. 1998c; Bah- 
cafl 1998; Bahcall et al. 1997; Fan ct al. 1997; Bartelmann 
et al. 1998; Lineweaver 1998), whether determined more 
or less independently or in combination with other param- 
eters, this prejudice is starting to weaken. Conservatively, 
these results can be summarised as 



P2{Xo,i^o) = L{no\OA,0.2). 



(16) 



A prior constraint on Qq is useful since lensing statistics 
alone, as expected and as our results show, cannot usefully 
constrain Oq- 

In addition, we also consider the product of P2(Ao, f^o) 
with the age constraint pi(Ao,rio)i 



P3(Ao,f2o) =pi(Ao,f2o) OP2(Ao,Oo). 



3. 3. General discussion of prior inforrn,ation 



(17) 



Using harsher constraints would mean that results would 
reflect almost exclusively the prior information as opposed 
to the information derived from lensing statistics. It is not 
the purpose of this paper to do a joint analysis of several 
cosmological tests, ^ but rather to examine lens statistics as 
a cosmological test. For practical reasons, an upper limit 
on f2o and upper and lower limits on Aq are required. On 
the other hand, it is sensible to combine the results with 
conservative constraints from other well-understood cos- 
mological tests where there is general agreement and little 
room for debate. Within our upper and lower limits, we 
present our results both with and without the constraints 



After this was found to conflict with too many observations, 
the prejudice against a non-negligible value of Ao weakened, 
and the new prejudice has been in favour of a flat universe 
with Ao -I- f2o = 1. 

but see Sect. 5 



discussed above. The density values and confidence con- 
tours of the three prior probability density functions are 
shown in the right column of Fig. 1. 

4. Calculations 

Following K96, we use the Hinshaw & Krauss (1987) soft- 
ened isothermal sphere model for modeling the light de- 
fiection properties of the lens galaxies. For this model, the 
lens equation reads 



x-y 



bx 



S + y/x^ + S^ 



(18) 



where x is the angular position in the lens plane, y the an- 
gular position in the source plane, b = 47r(cr/c)^ (-Dds/^^os), 
(J denotes the one-dimensional velocity dispersion of the 
dark matter, s denotes the core radius, s = s/Dod is the 
angular core radius and Dod, Dqs and £>ds denote the an- 
gular size distances between the observer and the lens 
galaxy, the observer and the source and the lens galaxy 
and the source, respectively. Still following K96, we model 
the distribution of elliptical and lenticular lens galaxies us- 
ing Schechter functions with constant comoving density 



rie = 0.61 ± 0.21 10"^ Mpc"^ 

{h = Ho 10~^ km~^ sMpc) and slope 



Ue = -1.0 ±0.15. 



(19) 



(20) 



The lens galaxy luminosities are converted to the dark 
matter velocity dispersions of the softened isothermal lens 
model by means of Faber-Jackson type relations. 



where 

7e = 4.0 ± 0.5 



and 



225.0 ±22.5 km s" 



(21) 



(22) 



(23) 



The core radii of the softened isothermal lens model are 
varied with the dark matter velocity dispersions according 

to 



(24) 



where e = 2.8 and s*e = 10/i~^pc. We consider ellip- 
tical and lenticular lens galaxies only. For the number- 
magnitude counts of quasars, we adopt the best-fit model 
from K96. We neglect here evolution, dust and other pos- 
sible systematic effects and refer the reader to K96 for a 
discussion. 

In our first calculations we apply Eq. (8) and compute 
the a priori likelihood 



piD\Xo,ilo,io) 



(25) 
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-4. -2. 2. 








-4. -2. 2. 

^) 

Fig. 1. Le/t column: The cosmological parameter plane. The four curved lines in the plot at the left are the isochrones 
toHo = 0.5, . . . , 0.8. The straight line marks spatially flat world models. In the white region, the antipodal redshift falls be- 
low z = 4.5, the redshift of the most redshifted multiply imaged object currently known (Gott et al. 1989; Park & Gott 1997). 
Right column: The prior probability distributions pi(Ao,f2o) (top panel), ■p2{\o,0,o) (middle panel) and ■p[i(\o,0,o) (bottom 
panel). The pixel grey level is directly proportional to the probability density ratio, darker pixels reflect higher ratios. The pixel 
size reflects the resolution of our numerical computations. The contours mark 0.61, 0.26, 0.14 and 0.036 of the peak likelihood 
for the parameters Ao and Qq, which would correspond to the boundaries of the minimum 0.68, 0.90, 0.95 and 0.99 confidence 
regions if the distribution were Gaussian 
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and the posterior probability density functions 
Pi(Ao,f^o|i?) =Pp|Ao,r!o,^o)®Pi(Ao,r!o), (26) 
P2{\o,no\D)=p{D\Xo,no,^o)<^P2{Xo,^o), (27) 
and 

P3(Ao,Oo|i?) =p(i?|Ao.fio.^o)$5'P3(Ao.fio) (28) 

in the limit where all nuisance parameters take precisely 
their mean values. To obtain an impression of the con- 
sequences of neglecting the uncerntainties of the nuisance 
parameters, in our second calculation we increase the value 
of the most uncertain nuisance parameter. Tie, by two stan- 
dard deviations. 

For the computation of the innermost integral on the 
right side of Eq. (3), we consider the detectability of im- 
ages in pairs: If the separation between the two closest 
images - these are always images 2 and 3, counting from 
the outside in - is more than the lower limit of the sur- 
vey resolution limit S{y), we define the image separation 
and flux ratio for the purpose of sample selection based on 
the two brightest images, usually 1 and 2. Otherwise we 
construct one image from the combined fluxes and flux- 
weighted positions of images 2 and 3 and define the image 
separation and flux ratio for the purpose of sample selec- 
tion based on this combination image and image 1. 

In general, if the separation between images 1 and 2 
is too large for the survey and the separation between 
images 2 and 3 is large enough, then the image separation 
and flux ratio for the purpose of sample selection should be 
based on images 2 and 3. However, the present surveys are 
sensitive to the largest separations due to isolated galaxies, 
so this case doesn't need to be addressed in this paper 
(i.e. implementing it would lead to the same results in the 
present case). 

For the calculation of the probabilities p{mi, Zi,9i) 
the function C{y) selects only those image configurations 
whose separation is ±10 per cent of the observed separa- 
tion 9i. 

Each of the three integrals on the right side of Eq. (3) 

is approximated to an accuracy better than 0.004 by a 
family of recursive monotone stable formulae (Favati et al. 
1991a,b). 

5. Results and discussion 

5.1. Information content 

Given some observational data D, some model parame- 
ters <p, and some prior and posterior probability density 
functions p{(p) and p{4>\D), the amount of information ob- 
tained from the data (e.g. Bernardo & Smith 1994) (on a 
logarithmic scale) is 

\p{ct>\Dy 



log[I{D)] = / ;.(,^)log 



P(0) 



dct>. 



(29) 



The amoLmts of information obtained from our sample 
data are given in the caption of Fig. 3. 



5.2. Results 

The left panel of Fig. 2 shows the constraints on the cos- 
mological parameters Aq and fio based only on the infor- 
mation obtained from the lens statistics. 

Quite good constraints can be placed on Aq, more or 
less independent of Qq- It is a well-known fact (see K96 
and references therein) that lensing statistics can provide 
a good upper limit on Aq . While in the past this has mainly 
been discussed in the context of flat cosmological models, 
it is of course more general (Carroll et al. 1992; Falco 
et al. 1998). Although no unexpected effects are seen, it 
is important to note that this is the first time Aq and Qq 
have been used as independent parameters in conjunction 
with a non-singular lens model in an analysis of this type. 

Our analysis shows for the first time that gravitational 
lensing statistics can place a quite firm lower limit on Aq as 
well, again more or less independent of Qq. The constraint 
is not as tight since the gradient in the probability den- 
sity is not as steep towards negative Aq as towards positive 
Aq. If this lower limit can be improved enough, it could 
provide an independent confirmation of the detection of 
a positive cosmological constant (see Sect. 5.3). On the 
other hand, this might be difficult, since Poisson errors in 
the number of lenses and uncertainties in the normalisa- 
tion of the luminosity density of galaxies introduce rela- 
tively large uncertainties in this region of parameter space 
(K96, Falco et al. 1998). The latter effect is illustrated in 
the right panel of Fig. 2, where rio, the galaxy luminosity 
density normalisation, is increased by two standard devi- 
ations: the derived lower limit on Ao changes much more 
than does the upper limit. Nevertheless, our robust lower 
limit is much better than the —7 mentioned in Carroll 
et al. (1992). 

Our results place no useful constraints on Hq. It is 
interesting to note the fact, however, that likely values 
of Ao and CIq are positively correlated. This is similar to 
most cosmological tests, a notable exception being con- 
straints derived from CMB anisotropies (see Sect. 5.3). 
Fortunately, constraints on CIq from other sources are quite 
good (Sect. 3.2). Often, this is cast in the form of a con- 
straint on Qo — Ao (e.g. Cooray et al. 1999) or, perhaps 
more practical, Aq — fio- This is a reasonable way or re- 
ducing the information to one number, at least when one 
is concerned with upper limits on Aq (or Aq — Oq) in a rel- 
atively low-density universe. Besides the obvious depen- 
dencies on confidence levels and assumptions made, when 
comparing constraints on Aq from different investigations 
one should keep in mind whether they are approximations, 
like Ao — flo in lensing statistics, and whether a value for a 
particular scenario (for example, for a flat universe) is the 
'obvious' definition or in fact describes the intersection 
of the k = line with the corresponding 2-dimensional 
confidence contour, which in general will give a different 
number. Also, some authors plot 'real' confidence con- 
tours while some actually plot contours at values which 



8 



R. Quast & P. Helbig: Gravitational lensing statistics with extragalactic surveys. I 



1.5 



0.5 




1.5 



° 1 



0.5 




-4. 



A) 



Fig. 2. Left panel: The likelihood function p(D|Ao, Ho, ^o)- AH nuisance parameters are assumed to take precisely their mean 
values. The pixel grey level is directly proportional to the likelihood ratio, darker pixels reflect higher ratios. The pixel size 
reflects the resolution of our numerical computations. The contours mark the boundaries of the minimum 0.68, 0.90, 0.95 and 
0.99 confidence regions for the parameters Ao and flo- Right panel: Exactly the same as the left panel, but the parameter rie is 
increased by two standard deviations 



would correspond to certain confidence contours were the 
likelihood distribution in the parameter space in question 
Gaussian. 

The left plot in the top row of Fig. 3 shows the joint 
likehhood of our lensing statistics analysis and that ob- 
tained by using conservative estimates for Hq and the age 
of the universe (see Sect. 3.2). Although neither method 
alone sets useful constraints on f2o, their combination 
does, since the constraint involving Hq and the age of the 
universe only allows large values of CIq for Xq values which 
are excluded by lens statistics. Even though the 68% con- 
tour still allows almost the entire Oq range, it is obvi- 
ous from the grey scale that much lower values of Qq arc 
favoured by the joint constraints. The upper limit on Aq 
changes only slightly while, as is to be expected, the lower 
limit becomes tightcir. Also, the change caused by increas- 
ing Tie by 2 standard deviations is less pronounced, with 
regard to both lower and upper limits on Aq, as demon- 
strated in the right plot in top row of Fig. 3. 

The middle row of Fig. 3 shows the effect of including 
our prior information on fio (see Sect. 3.2). As is to be 
expected, (for both values of Ue) lower values of fio are 
favoured. This has the side effect of weakening our lower 
limit on Aq (though only slightly affecting the upper limit). 

We believe that the left plot of the bottom row of Fig. 3 
represents very robust constraints in the Ao-flo plane. The 
upper limits on Aq come from gravitational lensing statis- 
tics, which, due to the extremely rapid increase in the op- 
tical depth for larger values of Ao, arc quite robust and 
relatively insensitive to uncertainties in the input data 
(compare the left and right columns of Fig. 3) as well 
as to the prior information used data (compare the upper, 
lower and middle rows of Fig. 3). The upper and lower 



limits on Qq are based on a number of different methods 
and appear to be quite robust, as discussed in Sect. 3.2. 
The combination of the relatively secure knowledge of Hq 
and the age of the universe combine with lens statistics to 
produce a good lower limit on Aq, although this is to some 
extent still subject to the caveats mentioned above. 

If one is interested in the allowed range of Aq, one can 
marginalise over to obtain a probability distribution 
for Aq. This is illustrated in Fig. 4 and Table 2. 

5.3. Comparison with other results 

For comparison with other results, as a first step one can 
examine the allowed range of Aq for the current 'best- 
fit' value for Qq, which we take, based on the work cited 
in Sect. 3.2, to be Uq — 0.3. (A more conservative esti- 
mate is reflected by using the prior probability distribu- 
tion P2{Xq, Oo) = L(Oo|0.4, 0.2) as shown by the dark grey 
curve in Fig. 4 and in Table 2.) On the other hand, previ- 
ous limits on Aq have often been quoted for a flat universe 
(K96 and references therein). We consider both cases in 
Tables 3 and 4. 

We do not do a comparison for the special case Aq = 

sinc;e this analysis of gravitational lensing statistics does 
not usefully constrain (any limits coming only from the 
prior information on Oq)- 

It is beyond the scope of this paper to do a full com- 
parison of different cosmological tests. Except for a few 
general comments, we therefore restrict ourselves to com- 
ments on the similarities and differences between the re- 
sults from this work without using prior information on 
Ao and Qq, i.e. (the left plot in) Fig. 2, and the those 
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Fig. 3. Left column: The posterior probability density functions p\{\o,Q,o\D) (top panel), p2(Ao, fiolf) (middle panel) and 
P3{\o,Q,o\D) (bottom panel). All nuisance parameters are assumed to take precisely their mean values. The pixel grey level 
is directly proportional to the likelihood ratio, darker pixels reflect higher ratios. The pixel size reflects the resolution of our 
numerical computations. The contours mark the boundaries of the minimum 0.68, 0.90, 0.95 and 0.99 confidence regions for 
the parameters Ao and Q.o. The respective amounts of information (Eq. (29)) obtained from our sample data are /i = 1.74, 
I2 = 1.24 and I3 = 1.74. Right column: Exactly the same as the left column, but the parameter Ue is increased by two standard 
deviations 
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Fig. 4. Left column: The top panel shows the normalised marginal likelihood function p{\o \D) (light grey curve) and the marginal 
posterior probability density functions pi(D|Ao) (medium grey curve), P2iD\\o) (dark grey curve) and p3(-D|Ao) (black curve). 
All nuisance parameters are assumed to take precisely their mean values. The bottom panel shows the respective cumulative 
distribution functions. These can be used to construct any desired f2o-averaged upper or lower limits on Ao- Right column: 
Exactly the same as the left column, but the parameter is increased by two standard deviations 



Table 2. Marginal mean values, standard deviations and 0.95 confidence intervals for the parameter Ao on the basis of the 
marginal distributions shown in the top row of Fig. 4; 'information' refers to Eq. 29 



Distribution 


Mean 


standard deviation 


95% c.l. range 




information 


p{D\\o) 


-0.35 


1.07 


-2.55 


1.51 




Pi(Ao|0) 


-0.02 


0.80 


-1.59 


1.50 


1.74 


P2{Xo\D) 


-0.78 


0.97 


-2.85 


0.76 


1.24 


P3{Xo\D) 


-0.34 


0.67 


-1.72 


0.79 


1.74 



from K96 and Falco et al. (1998) (using only optical data, 
i.e. the lower left plot in their Fig. 5). 

Taking all results at face value and examining the 
flo = 0.3 case first, we note that with 'three-and-one- 
half ' exceptions (counting as one test each the four from 
this work and the three from Falco et al. (1998)) the 68% 
c.l. lower limit from Lineweaver (1998) is higher than all 
68% upper limits from other tests, while the 95% lower 



and upper confidence levels from Lineweaver (1998) are 
higher than the corresponding limits from the other tests 
for all but one of these. Even at the 99.9% confidence 
level (not shown in Table 3), the Lineweaver (1998) re- 
sult requires Aq > 0.12. If one assumes f^o — 0.3, only 
Lineweaver (1998) requires Aq > 0, though all other tests 
(except Carlberg (1998)) are compatible with this. This is 
not surprising, since it is well-known that constraints from 
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Table 3. Mean values and ranges for assorted confidence levels for the parameter Ao for our a priori and various a posteriori 
likelihoods from this analysis and from other tests from the literature (using the latest publicly available results) for the special 
case Qq = 0.3. Except where noted, the ranges quoted are the projections of the corresponding confidence contours in the Ao-f2o 
plane onto the Ao axis" (as opposed to f2o-independent estimates, which of course would always give a smaller range), and are 
of course two-sided, not one-sided, bounds. Values are either those quoted in the references given and/or obtained from figures 
in those references; inequalities mean that the corresponding confidence contour is to be found in the range indicated by the 
inequality, e.g. < —1.2 would mean that the corresponding contour level is to be found at Ao < —1-2, not that the constraint 
is Ao < —1.2 at the corresponding confidence level. This arises because the corresponding area of parameter space was not 
examined in the reference in question. If the confidence interval could not be determined from the reference, both values in the 
corresponding column are missing 



Cosmological test 




68% c.l. 


range 


90% c. 


.1. range 


95% c.l. range 


99% c.l. 


range 


this work, p{D\\o) 




-1.18 


0.24 


-2.19 


0.50 


-2.81 


0.60 


-4.16 


0.73 


this work, pi(A()|D) 




-0.97 


0.46 


-1.55 


0.60 


-1.89 


0.69 


-2.73 


0.81 


this work, p2{Xi)\D) 




-2.00 


0.49 


-3.33 


0.65 


-4.10 


0.72 


< -5.00 


0.80 


this work, p:i(X[)\D) 




-1.20 


0.52 


-1.98 


0.69 


-2.35 


0.77 


-3.40 


0.86 


lens statistics (K96) 








not possible 


since only k 


= models 


considered 






radio lenses'"^ 




-0.54 


0.28 


< -1.00 


0.75 


< -1.00 


0.89 






optical lenses "^"^ 


< 


-1.00 


0.37 


< -1.00 


0.75 


< -1.00 


0.89 






radio -|- optical lenses-'^^ 


< 


-1.00 


-0.12 


< -1.00 


0.50 


< -1.00 


0.70 


< -1.00 


0.89 


supernovae m-z relation A'^ 




-0.70 


0.50 


-1.15 


0.75 










supernovae m-z relation jB'-'* 




0.78 


1.00 






0.53 


1.27 


0.27 


1.41 


CNOC survey' 


< 


-0.50 


< -0.50 


< -0.50 


< -0.50 










CMB™" 




0.44 


0.67 






0.36 


> 0.90 


0.26 


> 0.90 


CMB -t- IRAS° 








not possible 


since only k 


= models 


considered 






double radio sources'" 




0.00 


1.00 


< -2.00 


1.39 











" Note that some references quote confidence ranges for k = — in general, these will be different than the projection of the 
intersection of the corresponding contour in the Ao-Hq plane onto the Ao-axis. 
' Falco et al. (1998) 
" contour at 95.4% not 95% 

Falco ct al. (1998) 
" contour at 95.4% not 95% 
^ Falco et al. (1998) 
contour at 95.4% not 95% 

Perlmutter et al. (1998) 
' Riess ct al. (1998) 
^ Fig. 6, solid contours 

* contours at 68.3%, 95.4% and 99.7% instead of 68%, 95% and 99% respectively 
' Carlberg (1998) 
Lineweaver (1998) 

" contours at 68.3%, 95.4% and 99.7% instead of 68%, 95% and 99%, respectively 
" Webster et al. (1998) 
Guerra et al. (1998) 



CMB anisotropics tend to run more or less orthogonal in 
the Ao-ilo plane to those from most other tests (e.g. White 
1998; Eisenstein et al. 1998b; Tegmark et aL 1998a,b). 

Examining the k = case, it is interesting to note 
that the 68% (90%) confidence level lower limit on Aq 
from Carlberg (1998) is higher than all of the 68% (90%) 
c.l. upper limits from all other tests except Guerra et al. 
(1998). Otherwise, with 'one-and-one-half exceptions all 
tests are compatible even at the 68% confidence level. If 
one assumes k = 0, then the evidence for Aq > looks 
convincing: at the 68% confidence level, again with 'one- 
and-one-half exceptions, all tests indicate Ao > 0; even at 
90% the evidence is still quite good, if one keeps in mind 



that the gradient towards smaller values of Aq is generally 
not as steep as towards larger values. 

Again taken at face value, neither the k = case nor 
the Clo = 0.3 case are compatible with all tests, even at the 
w90% confidence level. It appears the simplest solution to 
achieve concordance would be to have ~ 0.2, which is 
within the error on Oq discussed in Sect. 3.2. For k = 
this would imply Aq = 0.8, which seems to be ruled out, 
thus ruling out the flat universe altogether. For a non-flat 
universe, reducing Oq would, due to the CMB constraint, 
require a higher value of Ao, and thus make the Ao = 
case more unlikely, ruling out this special case as well. 
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Table 4. Mean values and ranges for assorted confidence levels for the parameter Ao for our a priori and various a posteriori 
likelihoods from this analysis and from other tests from the literature (using the latest publicly available results) for the special 
case fc = 0. Otherwise the same as Table 3, in particular the references are not listed in the footnotes to this table. X denotes 

the fact that tli<>r<> is no intersix'tioii of tli<> confidence contour witli the k = line 



Cosmological test 


68% c.l. 


range 


90% c.l. 


range 


95% c.l. 


, range 


99% c.l. 


, range 


this work, p(D\X[)) 


-0.68 


0.51 


< -1.00 


0.57 


< -1.00 


0.62 


< -1.00 


0.70 


tliis work, pi(A()|D) 


-0.09 


0.56 


-0.38 


0.64 


-0.57 


0.68 


-1.04 


0.81 


this work, p2(Ao|-D) 


X 


X 


0.09 


0.69 


-0.03 


0.73 


-0.28 


0.92 


this work, p3(Ao|-D) 


0.47 


0.48 


0.18 


0.67 


0.07 


0.70 


-0.14 


0.84 


lens statistics" 










< 0.00 


0.66 






radio lenses' 


-0.47 


0.56 


< -1.00 


0.72 


< -1.00 


0.80 


< -1.00 


0.85 


optical lenses" 


< -1.00 


0.56 


< -1.00 


0.72 


< -1.00 


0.80 


< -1.00 


0.87 


radio + optical lenses'^ 


-0.87 


0.43 


< -1.00 


0.60 


< -1.00 


0.69 


< -1.00 


0.78 


supernovae m-z relation A 


0.20 


0.60 


-0.05 


0.75 










supernovae m-z relation B"-^ 


0.74 


0.83 






0.61 


0.92 


0.50 


1.00 


CNOC survey 


0.85 


0.95 


0.81 


0.98 










CMB» 


< 0.00 


0.60 


< 0.00 


< 0.00 


< 0.00 


< 0.00 


< 0.00 


< 0.00 


CMB + IRAS'" 


0.47 


0.71 














double radio sources 


0.35 


1.00 


0.70 


1.00 











" value for fe = 0, not projection 
contour at 95.4% not 95% 

" contour at 95.4% not 95% 
contour at 95.4% not 95% 
" Fig. 6, solid contours 

f contours at 68.3%, 95.4% and 99.7% instead of 68%, 95% and 99% respectively 

" contour at 68.3% instead of 68%; other contours, and part of the 68.3% contour, lie partially in the k = +1 area of parameter 
space which was not examined for technical reasons in Lineweaver (1998) 
value for fe = 0, not projection 



On balance, a cosmological model with Aq « 0.3 and 
Oo ~ 0.25 seems compatible with all known observational 
data (not just those discussed here) at a comfortable con- 
fidence level. 

For a 'likely' Oq value of 0.3 we have calculated the 

likelihood with the higher resolution AAq = 0.01. This is 
shown in Fig. 5. From these calculations one can extract 
confidence limits which, due to the higher resolution in 

Aq, are more acc\.irate. These are presented in Table 5 and 
should be compared to those for p{D\Xo) from Table 3. 

Again, a full discussion of joint constraints involving 
discussion of possible sources of error for each test, as well 
as comparing the full contours in the Ao-fio plane, is be- 
yond the scope of this paper. However, quick comparisons 
would be aided were the results of all tests available in 
an easy-to-process electronic form (see below); such quick 
consistency tests would enable one to spot areas of incon- 
sistency much more quickly. Also, it should be emphasised 
that the projections onto the Ao-axis of the intersection of 
a particular confidence contour with the flo = 0.3 or fc = 
axis are generally not the same as the corresponding con- 
fidence interval for the fig = 0.3 or fc = special cases. 

For a flat universe, our 95% confidence level upper 
limit on Aq-Oq, i.e. the value of Aq where this contour 

crosses the fc = line, is Aq < 0.62. This is essentially the 
same as the Aq < 0.66 of K96, as was to be expected con- 



sidering we used essentially the same data and methods. 
Interpreted cautiously, one might conclude from this that 
the singular isothermal sphere model is a good approxi- 
mation as far as determining the cosmological parameters 
from lens statistics is concerned, as was assumed in Falco 
et al. (1998). Our 99% confidence level upper limit on Aq 
is Ao < 0.70. This is quite a tight upper bound on Aq and 
appears to be quite robust. 

Perhaps more interesting is the comparison with (the 
results using only optical data in) Falco et al. (1998). Al- 
though a detailed comparison is complicated by the dif- 
ferent plotting scheme and reducing the entire contour 
(or indeed grey-scale) plot to a few numbers throws away 
information, it is obvious that the plots are broadly sim- 
ilar. Our 68% contour is, for fig ~ Ij roughly parallel to 
the Ao-axis at Aq « —1. This is just at the edge of the 
Falco et al. (1998) plot, and as they provide no grey-scale, 
it is difficult to compare the lower limits on Aq. Thus, 
while our main goal was to explore a 'large enough' region 
of parameter space, comparison in the areas where there 
is overlap shows consistency, which strengthens our faith 
in the conclusions pertaining to areas of parameter space 
where there is no overlap. 

Recently, it has become quite fashionable to discuss 
joint constraints derived from a variety of cosmological 
tests. This has grown from plotting the overlap of likeli- 
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Fig. 5. Left panel: The Ukelihood function as a function of Ao for CIq — 0.3 and with all nuisance parameters taking their default 
values. Right panel: The same but plotted cumulatively. See Table 5 

Table 5. Confidence ranges for Ao assuming flo = 0.3. Unlike the results presented in Table 3, these figures are for a specific 
value of r^o and not the values of intersection of particular contours with the ilo = 0.3 line in the Ao-fio plane. These are more 
appropriate if one is convinced that fio = 0.3 and have been calculated using ten times better resolution than the rest of our 
results presented in this work. See Fig. 5 



68% c.l. range 


90% c.l. range 


95% c.l. range 


99% c.l. range 


-1.27 0.27 


-2.26 0.51 


-2.87 0.60 


-4.10 0.72 



hood contours (often in a space spanned by parameters 
other than Ao and fio) (e.g. Ostriker & Steinhardt 1995; 
Turner 1996; Bagla et al. 1996; Krauss 1998; White 1998) 
to full-blown joint likelihood analyses, both detailed theo- 
retical investigations of what will be possible in the future 
(e.g. Tegmark et al. 1998a, b; Eisenstein et al. 1998a, b) and 
more restricted analyses using present data (e.g. Webster 
et al. 1998). While in some cases it is quick and easy to cal- 
culate the likelihood as a function of Aq and fio given the 
data, for example for tests using the m-z relation, in other 
cases such as the present one it is a major programming 
and computational effort to do so. To aid comparisons, all 
figures from this paper are available in the form of tables 
of numbers at 

http : //mult i vac . jb . man .ac.uk: 8000/ceres 
/data_f rom_papers/lower_limit . html 

and we urge our colleagues to follow our example. We ap- 
plaud the fact that most results are now presented in the 
Ao-fio plane, as opposed to using other parameters such as 
qo or ritot = Ao -I- ilo- A further aid in comparison would 
be a uniform choice of axes. We prefer to plot Oq on the 
y-axis and Aq on the x axis since up/down symmetry is 
less fundamental than left /right symmetry and this mir- 
rors the fact that Hq has the physical lower limit ^0 = 
whereas no corresponding upper or lower limits for Ao ex- 
ist. Square plots with the same range would further aid the 



comparison. Of course, if all data are publicly available, 
then they can be re-plotted to taste. 

6. Summary and conclusions 

We have re-analysed optical gravitational lens surveys 
from the literature, using the techniques described in Ko- 
chanek (1996), for the first time allowing both the cos- 
mological constant Aq and the density parameter Qq to 
be free parameters while also using a non-singular lens 
model. We confirm the well-known results that gravita- 
tional lensing statistics can provide a good upper limit 
on Ao but are relatively insensitive to Qq- Wc have pre- 
sented the new result of a robust lower limit on Aq , which 
is a substantial improvement on previously known robust 
lower limits. Coupled with relatively conservative prior in- 
formation about the Hubble constant Hq, the age of the 
universe and the well-established value of fig, one can re- 
duce the allowed parameter space in the Ao-f^o plane to a 
small, finite region, which is similar to the area allowed by 
joint constraints based on many other cosmological tests 
(see Fig. 3). 

Using lens statistics information alone, at 95% con- 
fidence, our lower and upper limits on Ao — fio are re- 
spectively —3.17 and 0.3. For a flat universe, this cor- 
responds to lower and upper limits on Ao of respectively 
— 1.09 and 0.65. Keeping in mind the difficulties of a quan- 
titative comparison, this is in good agreement with other 
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recent measurements of the cosmological constant. This 
value was calculated from Table 5 and assuming a degen- 
eracy in Ao — flo as in Cooray et al. (1999) and Cooray 
(1999). For comparison, from Tabic 4, the corresponding 
value for the upper limit on Aq is 0.62 and the value from 
K96 is 0.66.5 

For detailed comparison of cosmological tests, one 
needs to compare confidence contours — calculated in the 
same, preferably in the 'real', way — in the same parame- 
ter space. Of course, this makes it difficult to meaningfully 
reduce the results of a given cosmological test to one or 
even a few single numbers. Unless a cosmological test is de- 
veloped which can measure Aq independently of any other 
parameters, there is not much point in quoting unqualified 
'limits on Aq'. 

Presently tentative claims of the detection of a posi- 
tive cosmological constant, if true, would rank among the 
great discoveries of cosmology. Even though there are se- 
rious difficulties involved, it seems worthwhile to be able 
to confirm this result by improving the lower limit on 
Aq derived from gravitational lensing statistics. Targetting 
the two primary sources of uncertainty calls for improving 
our knowledge of the normalisation of the local luminos- 
ity density of galaxies as well as increasing the size of 
gravitational lens surveys. As far as the latter goes, the 
CLASS survey (Browne et al. 1998; Myers et al. 1999) 
looks the most promising at the moment. In a companion 
paper (Helbig et al. 1999), we have shown that compa- 
rable constraints to the ones presented in this work can 
be obtained from the JVAS gravitational lens survey; this 
gives us hope that the much larger CLASS survey will 
offer improvement in this area. 

Cosmological tests which set tight upper limits on flo 
imply, for a flat k = universe, a value of Aq which is 
ruled out by lensing statistics. For a non-flat universe, 
many tests are indicating Aq > 0, and at present a cosmo- 
logical model with Aq ~ 0.3 and Qq » 0.25 seems compat- 
ible with all known observational data, with neither a flat 
universe nor a universe without a positive cosmological 
constant being viable alternatives. The simplest case, the 
Einstein-de Sitter universe with Aq = and f2o = 1, both 
flat and without a cosmological constant, had been aban- 
doned long before the new observational data cited in this 
work came to light (sec, e.g., Ostriker & Steinhardt 1995, 
and references therein); this trend has continued, with the 
next-most-simple cases also no longer viable. For Aq and 
Qq, wc have in a sense reached the least simple case; it 
will be interesting to see if this trend continues with re- 
gard to the other cosmological parameters, in particular 
those which can be measured by the Planck Surveyor mis- 



^ The value from Cooray et al. (1999) and Cooray (1999) 
is 0.79, but it should be noted this value (the same in both 
papers) is based on different surveys, namely the Hubble Deep 
Field and CLASS, respecively.) 
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Fig. A.l. Likelihood that the non-lenses in our sample arc not 
lenses. The contour levels mark changes of a factor of ten in 
the probability, which is also indicated by the grey scale, darker 
values corresponding to higher values 

sion. Larger gravitational lens surveys such as CLASS will 
be a step in this direction. 
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Appendix A: Getting a feel for it 

The likelihood of a given cosmological model for a given 
set of observational data, calculated using Eq. (3), is the 
result of the complex interplay of many factors. While this 
is necessary for a detailed analysis, it perhaps obscures 
the fact that the likelihood is basically the product of two 
terms, the likelihood that the non-lenses in our sample are 
not lenses (see Fig. A.l) and the likelihood that the lenses 
in our sample (see Fig. A. 2) have the observed properties.® 
The latter in turn is the result of two basic effects: the 
dependency of the volume element dV/dz on Aq and flo 
(see Fig. A. 3) and the dependency on the lensing cross 
section on Aq and f2o (see Fig. A. 4). One can also use the 
probability that the non-lenses in our sample are not lenses 

® It is interesting to note that the measurement of Ao by Im 
et al. (1997) (who obtain Ao = 0.64lQ'2g for a flat universe and 
thus a lower limit) essentially corresponds to Fig. A. 2 (though 
with a different sample of lenses) . Since lensing is a rare phe- 
nomenon, small-number statistics arc a source of concern. The 
advantage of a well-defined gravitational lens survey, as op- 
posed to using a 'sample from the literature', is that the (much 
greater) number of non-lenses in the sample also makes a con- 
tribution. A comparison of Figs. A.l and A. 2 hints that not 
taking the non-lenses into account would tend to favour a high 
value of Ao, as indeed found by Im et al. (1997). 



R. Quast & P. Helbig: Gravitational lensing statistics with extragalactic surveys. I 



15 



2. 



1.5 



0.5 








-4. -2. 



Fig. A. 2. Likeliliood tliat tlie lenses in our sample have the 
properties they do. The contour levels mark changes of a factor 
of ten in the probability, which is also indicated by the grey 
scale, darker values corresponding to higher values 
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Fig. A. 3. The volume element dV/dz at the typical lens 
gala^xy redshift Zd = 0.7. The contours indicate the fraction 
0.1, 0.2, ... 0.6 of the volume element in the limiting case of 
the de Sitter model (Ao = 1, f2o = 0). This is also indicated by 
the grey scale, darker values corresponding to a larger volume. 
For smaller redshifts the contours are more vertical (and fur- 
ther apart), for larger redshifts more horizontal (cf. Fig. 3 of 
Tegmark et al. (1998b) but note their swapped axes) 

(illustrated in Fig. A.l) to calculate the expected number 
of lenses in our sample (see Fig. A. 5), although obviously 
just counting the number of lenses does not make use of as 
much of the available information as does using Eq. (3). 
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et al. (1999) and Cooray (1999) where the number of lenses as 
a function of Ao and ilo has been calculated for the Hubble 
Deep Field and for CLASS 
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